Mesoscopic Tunneling Magnetoresistance 
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We study spin-dependent transport through ferromagnet/normal-metal/ferromagnet double tunnel 
junctions in the mesoscopic Coulomb blockade regime. A general transport equation allows us to 
calculate the conductance in the absence or presence of spin-orbit interaction and for arbitrary 
orientation of the lead magnetizations. The tunneling magnetoresistance (TMR), defined at the 
Coulomb blockade conductance peaks, is calculated and its probability distribution presented. We 
show that mesoscopic fluctuations can lead to the optimal value of the TMR. 

PACS Numbers: 75.70.Pa, 73.23.Hk, 73.23-b, 73.40.Gk 



Spin-polarized electron tunneling has become a very 
active area of research, motivated by its potential role 
in new electronic devices Has well as by the possibility 
of observing novel effects In ferromagnet/ferro- 

magnet single tunnel junctions (F-F, separated by an 
insulating barrier), the subject has been investigated for 
some time ||, their main feature is that the resistance 
depends on the relative orientation of the magnetization 
in the ferromagnets. This is characterized by the tun- 
neling magnetoresistance (TMR), defined as the relative 
change of the resistance when the magnetizations of the 
ferromagnets rotate from being parallel to anti-parallel. 

A richer variety of effects are observed in ferromag- 
net/normal-metal/ferromagnet double tunnel junctions 
(F-N-F) as well as in completely ferromagnetic double 
junctions (F-F-F), due to the interplay between charge 
and spin. These systems consist of two ferromagnetic 
leads separated by a metallic grain. If the capacitance of 
the grain is such that the thermal energy ksT is smaller 
than the charging energy Eq, then the Coulomb block- 
ade (CB) of electron tunneling occurs. Specifically, the 
conductance shows pronounced peaks as a function of 
an external gate controlling the electrostatic potential on 
the grain. How this charging effect modifies the TMR has 
been the subject of recent experimental and theoret- 
ical @-|l|] work. In F-F-F systems, the TMR is enhanced 
in the Coulomb blockade regime by higher-order tunnel- 
ing processes such as co-tunncling while in F-N-F 
systems the TMR is reduced ]l0| . The effects of quantum 
interference within the Coulomb blockade has not been 
considered in these systems as of yet. 

On the other hand, the quantum Coulomb blockade 
has been extensively investigated in the absence of fer- 
romagnets |l^,[l3|]. Quantum interference in these struc- 
tures necessarily entails mesoscopic fluctuations as the 
detailed shape of the grains cannot be controlled. Such 
effects become important for sufficiently low temperature 
or, at room temperature, for sufficiently small systems. 
In such a quantum regime (fcsT<C A), the CB conduc- 
tance peak height, G pea k, fluctuates strongly as a func- 
tion of the gate voltage 



two fields, quantum mesoscopics and spintronics. Specif- 
ically, we study quantum effects in the TMR at the CB 
peaks — the mesoscopic tunneling magnetoresistance — in 
F-N-F systems. In the case of zero spin-orbit tunneling, 
we show that there are CB peaks with the maximum pos- 
sible value of TMR; in contrast to the classical case, no 
special tuning is required. For strong spin-orbit coupling, 
the TMR can be large despite the SO coupling. As for 
Gp Ca k, the mesoscopic TMR is characterized by a prob- 
ability distribution rather than by a single number, and 
we obtain this distribution in both cases above. 

Since G pea k depends on the properties of a single wave 
function, one can legitimately wonder if it is possible to 
probe its spin structure (a very active subject ]5|, p^ , p^| ) 
through the conductance in non-collinear configurations. 
For such a configuration, we calculate G pea k as a function 
of the angle between the magnetizations in the two leads. 
In this case as well as for strong SO, the usual descrip- 
tion in terms of a rate-equation is not possible, and a 
more general equation for the conductance is presented. 
Finally, we study the case of broken spin degeneracy in 
the grain and point out that it is equivalent to having a 
half-metallic jll| lead. 

Figure 1 shows a schematic of the system: a nor- 
mal metal grain is weakly coupled to two ferromag- 
netic reservoirs by tunnel junctions. The potential of 
the grain is controlled by gate voltage V g . We consider 
r<fc s T<A<£;c, where T is the total width of the 
resonant levels in the dot. In this regime, only a sin- 
gle energy level contributes to the conductance. Because 
of time-reversal symmetry, each level is doubly degener- 
ate (Kramers doublet). Within a Hartree-Fock treatment 
of the electron-electron interactions, the transport is de- 
scribed in terms of the self-consistent single electron wave 
functions; for the resonant doublet, these are the spinors 
*i(r) = (0(r), X (r)) T and * a (r) = (-** (r), 4>* (r)) T . The 
mesoscopic tunneling magnetoresistance is defined as 
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Here we study a problem at the intersection of these 



where Gp^p) is the conductance at the CB peak in the 
parallel (antiparallcl) configuration. 
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FIG. 1. Schematic picture of the F-N-F double junction. 
A normal metal grain is coupled to ferromagnetic leads by 
tunnel junctions and is small enough for charging to be im- 
portant. When the magnetization of one lead is flipped, the 
conductance changes: the relative change with respect to the 
parallel configuration defines the TMR. 



Absence of Spin- Orbit — Consider the simplest case of 
no SO interaction and collinear magnetizations. Then, 
the Kramers doublet corresponds to two-fold spin degen- 
eracy, x( r ) = Oj assuming the Zeeman energy associated 
with any magnetic field present is much smaller than 
k R T. Linear-response theory yields for the conductance 
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where the sum is over the spin degenerate states and 
A is a numerical factor [fl7|— p~9|l . The partial width T^ n , 
q = L(R), due to tunneling to the left (right) electrode 
is = 2tt T, a P q J V am\ 2 with V nm the matrix element 
between \m) in the grain and the channels a=],{ in lead 
q, and p q a the local density of states of those channels. 
For point-contact leads, rL> = |0 g | 2 tL^ where |</> g | 2 is 
the wave function of the resonant level at the contact 
point r q and t q a is a factor which depends on both p q a 
and the tunnel barrier. Notice we assume that the tun- 
neling through the barriers conserves spin. Letting D q 
(d q ) denote t q a for the majority (minority) spins, we find 
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Since metallic grains are generally irregular in shape, 
one expects the chaotic nature of the single particle clas- 
sical dynamics to produce random matrix statistics for 
the wave functions. Following Jl2], we use RMT to de- 



scribe the fluctuations of I 



the Gaussian Orthogonal 



Ensemble (GOE) for the Hamiltonian implies the Porter- 
Thomas distribution, P(|0 ? | 2 ) oc exp(— |</> g | 2 ). As- 
suming no correlation between \4>l\ 2 and \(j> R \ 2 because 
the contact points are far apart, we obtain 
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a = (c 1/4 + r 1/4 )/(^ 1/4 + ^ 1/4 ), 

C = D R d R /D L d L , S = D R d L /D L d R , (5) 
Vm = 2P L P R /[1 - P L P R + [(1 - P 2 )(l - P 2 )] 1 / 2 ] . 



Here 



P, 
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L, R 
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is the spin-polarization of the electrons coming from the 
leads pc|] . Note that in the case of symmetric barri- 
ers and the same ferromagnets — and only in this case — 
Pgoe(??) factorizes into a term depending on the polar- 
izations P q and one depending on the properties of the 
grain. 

The maximum value of the TMR, r) m , is also the upper 
limit for the TMR when k R T » Ec- In that regime, it 
is reached when Q — 1 (optimal asymmetry) while in our 
case the condition is \4*l\ 2 /\<t>R\ 2 = Q- If the two leads are 
made of different materials, the former condition requires 
that the tunneling barriers be expressly designed while in 
the mesoscopic regime there are always fluctuations that 
satisfy the latter condition. Indeed, note that the prob- 
ability of the TMR being close to the maximum value is 
not small. In this case, then, the mesoscopic fluctuations 
lead to the optimal value of the TMR. 

Broken Spin Degeneracy — The spin degeneracy can be 
removed by applying a magnetic field such that the Zee- 
man energy exceeds k B T . Then, only one channel con- 
tributes to (||): the resonant level acts as an ideal spin 
filter rendering the polarization of the left lead irrelevant. 
In this way an effective half-metallic Jl6| injector, an ob- 
ject of considerable interest in the spintronics community, 
can be constructed in conventional materials. Proceeding 
as before, we find for the TMR 
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with 77 m = 
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Results for this case can be obtained from the spin degen- 
erate results by formally taking d^ — > 0. This yields the 
distribution Eq. (|J) and the mean value (77) =r) m /(l + a), 
as before. Note that for D R <§; Dl , the average TMR is 
much bigger than the symmetric single tunnel-junction 
result ||: (77) — rj m compared to 2P R /{\ — Pj|), a fac- 
tor (I + 1/Pr) larger. For ferromagnetic leads made of 
Co (P~0.35) this means an enhancement of 300%. The 
magnetic field here must, of course, be smaller than the 
coercive field of the electrode (in the AP configuration) , 
making low temperatures advantageous in observing this 
effect (typically, B < 500 G and T<70mK). 

Spin-Orbit Coupling — SO coupling causes the direc- 
tion of the spin of an electron to rotate while in the 
grain, so that the natural axis of quantization on the left 
is not the same as on the right. Time-reversal symme- 
try still applies but now neither component of the spinor 



2 



(</>(r), x(r)) T vanishes. Since the quantization axes in the 
leads are fixed by the bulk magnetization of the ferromag- 
nets, the electrons tunnel into a superposition of "up" 
and "down" states. This coherence cannot be included 
in a simple rate equation — a more general approach is 
required. Technically, the mixing of the spin channels 
means that T L and T R , which in the general case are 
matrices T^ m = 2irJ2 a Pa^?«K?m> cannot be diagonal- 
ized in the same basis. In the relevant 2x2 subspace T q 
takes the form 
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Note that [T L , T R ~\ ^ only when both leads are ferro- 
magnetic and SO coupling is present. 

From the Keldysh formalism, the current through an 
arbitrary system connected to two leads is [f2l|-p3"[ 



J = ^ Jds Tr { (T L - T R ) iG< + (T L f L T R f R )A} 

(9) 

where A = i(G r — G°) is the spectral function, G r ( a ) 
is the retarded (advanced) many-body Green function of 
the full system, and Jl(r) is the Fermi distribution in the 
left (right) lead. It is convenient to express G < , the lesser 
Green function, in terms of the operators I 9 = r 9 (G < — 

i/,A): g< - (r L + r*)- 1 ^ + i fi + i(r L /L + r R f R )A). 

For elastic transport and within the HF approximation, 
(I L + I R ) nn = 0, so that Eq. (§) becomes 

J = | fdeTr {T R (T l + r^) 1 T l a} (f L - f R ) . (10) 

This equation simplifies for linear-response and weak cou- 
pling. A is then evaluated in equilibrium: it depends 
on the equilibrium populations of the grain eigenstates 
(which depend only on the energy) and on the eigenfunc- 
tions of the isolated grain Q. Therefore, it can be shown 
that within the resonant subspace A is proportional to 
the identity. For the conductance, we finally get 
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Tr {r R (r L + t r ) 1 t l } (11) 



where the numerical factor A' has the same origin as A in 
(||). Eq. ( |lf"| ) contains (||) as a special limit; as far as we 
know, it is presented in this form for the first time. 

Since G pca k can be calculated in any basis, we choose 
one where xl — 0. Then, for symmetric leads and barriers 
(P q = P, C = l) the TMR from Eqs. (§) and (0) is 
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1 + (1-P 2 ) £-P 2 cos( 
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FIG. 2. Probability density of TMR for symmetric bar- 
riers and same ferromagnets with P — 0.35, 0.5, 0.8 in 
the symplectic ensemble (GSE). The TMR is scaled with 
?7 m = P 2 /(l - P 2 ). The distribution for P = 0.35 is multi- 
plied by 2 for purposes of comparison. The functional form 
of the distribution depends strongly on the value of P. In 
all cases, there is a cusp at n = 0. Note that rj can be both 
positive and negative and significantly different from zero. 
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with \^ q \ = |*i(r g )|. Note that -q = if the spin at the 
right contact is perpendicular to the polarization in the 
leads. For simplicity we will assume that SO is strong and 
so describe the statistics by the Gaussian Symplectic En- 
semble (GSE). The intermediate case is straightforward 
numerically. Then, </> and \ are uncorrelated complex 
random variables. Using their Gaussian distributions we 
find V{£) and show that "P(£, cos 9 S ) =V{£) /2. Therefore, 
the distribution of the TMR is 



0<r?<r? n 



(14) 



where 8 g is the azimuthal angle of the spinor ^(r^) and 



where A(x) = y/rj m — x (6 + 2?y m + x) and r] m = P 2 / (1 — 
P 2 ). The upper and lower limits correspond to the elec- 
tron tunneling out either conserving its original spin (r] m ) 
or with opposite spin (— P 2 ). Figure 2 shows the distri- 
bution for three different values of P. In the asymmetric 
case, the shape of the distribution is strongly dependent 
on the ratio D^/D R (not shown). 

We wish to emphasize three aspects of Vgse(v) : ( a ) 
The TMR can be either positive or negative for any po- 
larization as a consequence of the spin-rotation inside 
the grain— note that 7 5 gse('7>0)=7'gse(?7<0). (b) The 
functional form is strongly affected by SO — note in par- 
ticular the disappearance of the divergence at the maxi- 
mum value and the cusp at zero (which only vanishes for 
P— > 1). (c) Even for strong SO the TMR can be close 
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to its maximum zero-SO value] this is completly differ- 
ent from the situation in the classical regime, where the 
TMR in the presence of SO is always smaller than in its 
absence fl(| (note that the GSE represents the strongest 
SO coupling). This robustness of the mesoscopic TMR 
against SO is a consequence of tunneling through a single 
level. 

Non-Collinear Magnetizations — If the magnetizations 
of the two ferromagnetic leads are not collinear, a rate 
equation does not apply (even in the absence of SO) 
for the same reason as in the SO case — tunneling in- 
volves a coherent superposition of spin states. Instead, 
we must use Eq. (|ll|). We assume the magnetization on 
the left lead points in the z-direction while on the right 
lead it is tilted by an angle 29. The spin channels on 
the right lead are then given by \n, t) = (cos6>,sin6') T , 
|n, J.) = (— sin#, cos#) T . Writing \&i(i\r) and ^(rfl) in 
this basis, e.g. ^(r^) = (cj)RCOs8 + xr sin 9) \ri, t) + 
(— (j)RsinO + xr cos @) |^>4)> we readily obtain Tr and 
then the conductance. Using the basis where xl = we 
get 



G'p ea k(20) — Gp C ak(2#m) 



1 + (1 - P 2 )£ - P 2 B(29 m ) 
1 + (1 - P 2 )£ - P 2 B{29) 



(15) 



with B{x) = cos x cos 9 g + cos ip g sin x sin 9 g . Here, 9 g and 
tp g describe the rotation of the spin inside the grain in go- 
ing from left to right. Notice that the maximum is shifted 
in the presence of SO — the maximum (or minimum) 
conductance occurs for tan(2# m ) = tan(6> g ) cos(^ g ) — and 
that it becomes narrowed as P— >1. It is worth mention- 
ing that Eq. (|l^) is different from the result obtained 
in the classical regime |24|. Fig. 3 shows G pG ak(2#) for 
different polarizations both in the absence and in the 
presence of SO. 

This angular dependence might be used to study the 
spin-structure of single wave functions. A spin-polarized 
STM [jl5| could provide a spin-resolved image of a single 
wave function on the surface of a small grain, carbon 
nanotube or quantum corral. In principle it might also 
be used for imaging ferromagnetic quantum dots. This 
potential for revealing the spin structure of a single wave 
function deserves further investigation. 
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